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Abstract 

We present here a survey of recent results concerning the mathematical analysis of instabilities 
of the interface between two incompressible, non viscous, fluids of constant density and vorticity 
concentrated on the interface. This configuration includes the so-called Kelvin-Helmholtz (the two 
densities are equal), Rayleigh- Taylor (two different, nonzero, densities) and the water waves (one of 
Qh ' the densities is zero) problems. After a brief review of results concerning strong and weak solutions of 

•^C the Euler equation, we derive interface equations (such as the Birkhoff-Rott equation) that describe 

the motion of the interface. A linear analysis allows us to exhibit the main features of these equations 
(such as ellipticity properties); the consequences for the full, non linear, equations are then described. 
In particular, the solutions of the Kelvin-Helmholtz and Rayleigh- Taylor problems are necessarily 
analytic if they are above a certain threshold of regularity (a consequence is the illposedness of the 
initial value problem in a non analytic framework). We also say a few words on the phenomena 
that may occur below this regularity threshold. Finally, special attention is given to the water waves 
problem, which is much more stable than the Kelvin-Helmholtz and Rayleigh- Taylor configurations. 
Most of the results presented here are in 2d (the interface has dimension one), but we give a brief 
description of similarities and differences in the 3d case. 
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1 Introduction 

This contribution is a survey of recent results concerning the mathematical analysis of instabilities in 
interfaces between two fluids of different densities and vorticity concentrated on the interface. In 2d this 
interface is a time dependent curve and in 3d a time dependent surface. This includes the two following 
extreme cases: 

• The density of one of these fluids is infinitely small with respect to the density of the other and 
therefore can be taken equal to 0; this corresponds to the study of water waves or droplets in 
vacuum. 

• The two fluids have the same density (in fact there is only one fluid) and the vorticity is concentrated 
on an interface. 

Experiments and numerical simulation show that the range of possible behaviors for the interface is 
very rich. In some cases (especially in the water waves case) the interface is "stable" and its description 
has led to many studies in both the physical and mathematical literature. In other cases (this is the 
generic case when both densities are nonzero), the interface is "unstable". Our purpose is to show how 
mathematical analysis can contribute to the understanding of such instabilities. To do so, one starts 
from the most abstract model which is the initial value (or Cauchy) problem for the incompressible Euler 
equation. It models a configuration where viscosity and surface tension are neglected and therefore it 
is already too abstract to describe physical situations. As recalled below it admits solutions that would 
correspond to instantaneous creation or extinction of energy and as such would solve the energy crisis. 
However since in many cases the surface tension and viscosity are small, it turns out to be essential to 
study this limit case. We will consider the evolution of a fluid in a domain f2 which can be either the 
whole space SI = M. d a sub-domain ft C M. d or a periodic box f2 = (R/Z) d . The natural quantities are 



the density p the velocity u(x,t) and the pressure p(x,t). The equation of motion in ft are then given 
by the Euler equations (or conservation of momentum) , 

dt(pu) 4- V- (pu®u) + Vp = p§, (1) 

where the gravity g 6 K d is a constant vector, and the conservation of mass, 

d t p + V • (pu) = 0. (2) 

The fluid is also assumed to be incompressible; together with the conservation of mass, this yields the 
condition 

V-u = 0, (3) 

which can be viewed as a constraint on the evolution of the fluid with the pressure as the Lagrange 
multiplier of this constraint. In the presence of boundary (SI ^ M. d ) one has for all times t the following 
impermeability condition , 

Va; G dtt u(x, t) ■ ft — with ft the exterior normal on dQ, boundary of fi. (4) 

If $7 is unbounded, it is also assumed that u(t, ■) vanishes at infinity. 

Besides the time variable t, two coordinate systems appear naturally. In the Eulerian system, u(x,t) 
describes the velocity field at the point x and at the time t. In the Lagrangian system, we rather use the 
quantity u{x, t) that stands for the velocity of the fluid particle that was at the point x at the initial time 
t = 0. Equivalently, u{x,t) = u(X{x 1 t),t) 1 where 1 1-> X{x,t) solves the ordinary differential equation: 

X(x,t) = u(X(x,t),t), X(x,t = 0)=x. (5) 

We consider solutions u(x,t) which are continuous away from an interface £(£) which defines a partition 
of 57 = S7±(i) U S(i). We assume that on both sides of this interface the density of fluid p — p± is 
constant. It is important to keep in mind that the velocity of the fluid is well defined on Cl±(t) and 
therefore that the equation ([SJ , together with the incompressibility condition ([3]) , defines two families of 
volume preserving diffeomorphisms from Q,±(t = 0) onto Q±(t). However the velocity has no intrinsic 
definition on £(£) . Only the global evolution of this interface matters; as a consequence there is some 
freedom in the choice of the parametrization r(t, A) of the interface £(£), and this choice will play an 
important role below. 




Observe now that as long as the solution is smooth, the vorticity w = V A u is solution of the system 
of equations: 

dtio + u ■ Vcj = u ■ Vu , V • u = , V A u = u, in Q. (6) 



In 2d with a planar flow, the vorticity can be written u — (0,0,diU2 — d2U\) T and the term u ■ Vu 
vanishes. Still denoting by w the scalar vorticity w = d\U2 — c>2iii, the vorticity equation becomes 

d t uj + u ■ Vcj = 0. (7) 

In any case, one has 

u(x, 0) = for x £ S(0) =*> u(x, t) = for a; £ S(t), (8) 

as long as the dynamics is well defined. Several contributions have been devoted to the case where w(0) 
is not away from S(0) (cf. [3^1 EH HH1 [53 [Ml)- However since the evolution of the interface is the 
dominant effect we will consider only dynamics with vorticity away from the interface. Therefore the 
only parameters are the Atwood number: 

a= P+~P- -1 < a < 1, (9) 

P++P- 

and the gravity g. The Rayleigh- Taylor problem corresponds to < \a\ < 1 <^> p + / p_ , p+p- ^ 0, the 
Kelvin Helmholtz to a = ■£> p_ = p+. Finally the situation where a = ±1 <^=> p_p + =0 corresponds to 
a situation where one of the two densities vanishes and it describes the dynamics of water waves. The 
Rayleigh- Taylor and the Kelvin-Helmholtz problems are highly unstable while the stability of the water 
waves problem is ensured if an extra condition {called in the literature Taylor Hypothesis or Rayleigh- 
Taylor sign condition, see 37. 3j) . The above names indicate that the study of these problems goes back 
to the end of the nineteenth century or the very beginning of the twentieth century. 
The following facts should be kept in mind: 

• The range of applications of the analysis is huge: From classical water waves (and application in 
oceanography) to plasma physic in laser confinement and astrophysics. 

• The basic properties of the instabilities have already been established by the forerunners of the 
theory (which carries their names) by modal analysis of linearized problem near stationary solutions. 
A systematic use of the modal analysis can be found in the book of Chandrasekhar [12] . 

• The more recent contributions aim at a more systematic theory and allow for solutions that are 
not in a small neighborhood of stationary solutions. 

• We are interested here in the appearance of singularities and focus therefore on the most singular 
component of the solution (this is the central idea of symbolic analysis). However, the most singular 
parts are not necessarily the most important for the applications (in the so called shallow water 
regime the terms neglected by the symbolic analysis are the most important ones, see §7.5| . 

For the sake of simplicity most of the presentation is devoted to 2d problems. Similarities and differences 
which appear with the 3d configurations are considered in Section [5] 

2 Strong and weak solutions of the 2d Euler equation 

When the density of both fluids is the same (say, p+ = p- = 1, and thus Atwood number a = 0), the 
problem is reduced to the "genuine " incompressible Euler equations dlj) , (JSJ) . In this case, the gravity 
does not play any significant role (it can be put into the pressure term), and we thus take g = in (fT]): 

d t u + V • (u (8> u) + Vp = 0, V-u = 0. (10) 

By standard methods [40] one can show the following local existence theorem: 



Theorem 1. For any initial data uq G C 1,Q (R d ) (a > 0, d — 2, 3) with V • uo = on M. d , there exists a 
finite time T* > such that, with this initial data, one has a a unique classical solution to \ 10)) on the 
time interval [0, T*] . 

The above result is only local in time. T* depends on the size of the initial data and to the best 
of our knowledge this is the only "general" result. Since we are concerned with solutions having their 
vorticity concentrated on a curve £(£) (in 2d, or a surface in 3d), the above theorem does not help much 
and weak solutions have to be considered (for a precise meaning of what is meant by "weak solution" , 
see for instance [40]). Using the conservation of the vorticity along Lagrangian trajectories implied by 
(O, one gets the following classical results: 

Theorem 2. 1- (Yudovich 1621) With an initial vorticity uj G L°°(R 2 ) n L 1 (K 2 ), the Euler equations 
HU\) have a unique global weak solution which depends continuously on the data. 

2- If the initial vorticity is a bounded signed measure fl8f or with a finite number of changes of sign \39^ , 
the Euler equations hlU\) have at least one weak solution u(x,t) G C wea fc(Rt,-£ 2 (R 2 ))- 

However one should observe the big "gap" between the hypothesis of the point 1 and of the point 2 (this 
latter includes configurations allowing concentration of vorticity on an interface) . Without the hypothesis 
us G L°°(IR 2 ) no uniqueness or stability result is available. Furthermore wild solutions, for instance with 
space time compact support have been constructed by Scheffer [48] . Shnirelman [55] and finally by De 
Lellis and Szekelyhidi [T7] (see also the nice review paper by Villani [SB]) who have obtained the following 
instability, non uniqueness theorem (note that the second point implies that the wild solutions are not 
that exotic, since they are obtained as limits of strong, smooth solutions of the Euler equations): 

Theorem 3. 1- For any bounded domain O C Rj x R t (d = 2, 3) there exists a weak solution u G 
L°°(W^ xl f ) to flty) which satisfies the relation 

(x,t) eO^ \u{x,t)\ = l,(x,t)<£0^ \u(x,t)\ = 0. 

2- Moreover, u is the limit in L 2 (M. d x R t ) of a sequence of smooth, compactly supported, solutions u^ 
(k G N) to the Euler equations with smooth (and vanishing as k — > ooj source terms fk : 

d t u k + V • (u fc <g> u k ) + Vp fe = / fc , V • u k = 0, f k ^0 in V'(R d ). 

3 Weak solutions and interface equation 

In this section we derive a closed system of equations for the interface £(£) when the vorticity is a density 
concentrated on this curve, assumed that it satisfies the following smoothness assumption: 

Definition 4. 1- A rectifiable curve £ parametrized by its arc length, £ = {r(s)}, is chord-arc if there 
exists two constants < m < M < oo such that one has: 

m\si - s 2 \< \r(si) - r(s 2 )\ < M\s\ - s 2 \. 

A local version (near any point of Ti) of this definition is obtained with minor modifications. 

2- We say that a curve £ C £! is regular if it is bounded away from dCl, rectifiable and chord-arc . 



3.1 Biot-Savart formula and Hilbert Transform 

From the incompressibility relation V • u = one infers the existence of a stream function if) such that 

u(x,t) = V ± ip{x,t) = (-dyip{x,t),d x il>{x,t)) T . (11) 

Classical results of potential theory show that u(x, t) converges to two values u±(x, t) on the two sides of 
the interface. The incompressibility condition implies the continuity across the interface of the normal 
component of the velocity, but not of the tangential ones. Denoting by u?j_ the trace of the tangential 
component of the velocity on both sides of the interface, and by r(-, t) a parametrizatiorQ of the interface, 
£(£) = {r(A, £)}, standard distribution calculus shows that the scalar vorticity w is a density measure 
supported on S(t), 



V^eC(ft), (<p,w)=l <p(r(X,t))Z(X,t)d\, with u(X,t) = -\d x r\(u T + -uL) hr ^ tht) ; 

the quantity ui is called the vorticity density associated to the parametrization r(A,£). The vorticity 
density associated to the arc length parametrization is called the vortex strength, and will be denoted 
"f(s,t) throughout this text, 

7(g(A, *),*)= ,- /. ,xx. S(A,t), with s(X,t)= f \d x r(X',t)\dX'. (12) 

\d\r{X,t))\ J 

In terms of rp, one has 

A^ = w = w(8)5 E ( t ). (13) 

From (fTTj) . (fl3|) and the fact that u vanishes at inhnity if D, is unbounded, one can express u(x,t) for 
all x £ S(4) in terms of the vorticity density w by a Biot-Savart type formula. Denoting by G(x,y) the 
Green function of the Laplacian in fl with Dirichlet boundary condition, one has 

V* G fl(t)\E(t), «(»,<) - / V±G(x,y)Lj{t,y)dy 

V^G{x,r(X,t))u{X,t)dX. (14) 

For a; £ f2\S(i) the function u(t,x) is very smooth (indeed analytic). Though this is no more true for 
X 6 £(£), the formula (fT4l) can be extended to the case x £ S(4) provided the left hand side is replaced 

by 

u + + u_ 
w = 



2 ' 
and that on the right hand side the integral is defined as a principal value: 



v(x,t) =p.v. / V ± G(x,r(A,t))w(A,i)dA. (15) 



1 It is of course assumed that the interface £(t) is an orientable curve. We choose the orientation such that the normal 
vector n points from S7_(t) to 0+(t). The orientation of the tangent vector r is chosen such that (r, n) is a direct 
orthonormal basis of IR 2 

2 Changing A into —A if necessary, we always assume that d\r is positively colinear to r, that is, f = rs — rd\r. 



In particular for O = M 2 one recovers the standard Biot-Savart formula: 

v{l , t) ,l-R, v ,.j±Zl^ a( X, ()d X, ,16) 

with Re. denoting the rotation of angle ^ . Furthermore, due to the ellipticity of the Laplacian and also 
for fi ^ R 2 , a local version of the above formula remains valid modulo a more regular term 

v(x,t) = ^-%p.v. / X ~ r ^ u(X,t)dX + Ru(wy, (17) 

2tt 2 J{ r (x,t)eu} \x-r{X,t)\ 2 

in (fTTj) . a; i-> Ru(uj)(x) is an analytic function and w >->• Ru{oj) is a "regularizing operator " which 
depends on f2 and the open set 14. 

Remark 1. Denoting by v T and o n the tangential and normal components of v, one checks that: 



1 f „ z-r(A,i) 

— P-v. / n ■ -p 

2tt J{ r (\,t)eU} k — r(A,t)| J 



v T (x, t) = -— p.v. / « • ,- „,\'' 5(A, i)dA + i£( w ) (18) 



t,"(M) = ^p.v. / f • * ^ 5(A,f)rfA + J%M. (19) 

By a Taylor expansion of r(-,i) in the neighborhood of x, one checks that (x— r(X,t))-n = 0(|#— r(A,£)| 2 ) 
while one only has (x — r(X,t)) ■ r = 0(\x — r(X,t)\). It follows that v T is "one degree " more regular 
that v n . 

One of the main issues in the analysis of the Ray leigh- Taylor (and Kelvin-Hclmholtz) problem is the 
ellipticity of the linearized operator. This property is local by its definition and its effects. Therefore 
it is convenient to write, with a Galilean transform, the equation (|17|) in the neighborhood of a point 
r(A, t) = (0, 0)), with r, n given by the coordinate axis, and where £(£) is a graph (xi, ecr(xi, t)) which is 
a small perturbation of the x\ axis in a neighborhood \x\\ < S of the origin. In this setting (|17[) becomes 



Vl (x,t) = -^p.v. / 7 "7' 2 PiP - ( , m Co(x' x ,t)dx' x , (20) 

V 2 {x,t) = ^-p.v, I — 3737^ _ -7-7 ,^ 2 ^1,^1, (21) 



_e_ / e{xi,t) -a(x[,t) 

2ir P ' V ' J\ Xll<s ( Xl - x[) 2 + e 2 (a{x!,t) - a(x[,t)) 
If X\ — x'i 

2^ ' ' J\ X1 \<8 (xi -x[) 2 + e 2 (tr(xi,t) -a(x[,t)) 2 

where we omitted the more regular terms coming from Ru(uj). We thus see that two operators will play 
an important role in the analysis of the linearized equation (see []4] below): the Hilbert transform H and 
\D\, defined respectively as 

Hf{ Xl ) = -p.v. /^— r /( a; ' 1 )d a: ' 1 = J F- 1 (_i8 g n(0/(0), (22) 

77 J X\-X x 

mix,) = -p.v. / /M^MJ =dxi{Hf{xx)) = F -i m f i0) , (23) 

3.2 Jump relations and equations on the interface 

With minimal regularity hypothesis on the interface and on the vorticity density one can deduce "jump 
relations" . We have already seen that incompressibility relation V-u = and the identity VAu = 2®<5wt) 



lead to 

[«»] = and [ u t] = -±-Z, (24) 

where we use the standard notation [u] = u + — v,-. 

From the conservation equations of mass @ and momentum ([1}, written in the sense of distributions, 
we deduce further two equations for the evolution of the interface. If r(-,t) is a parametrization of the 
interface, £(£) = {r(X,t)}, one gets (cf. [54], and [14] for a slightly different, but equivalent formulation): 

(d t r -v)-n = (25) 

d t (-uj - a\d\r\v T ) + d\(-^ — - (-2 - a\d\r\v T )(v - d t r) ■ fj 

Observe that the above system does not fully determines the evolution of the vector r(X, t). The equation 
([25| concerns only the normal component of its time derivative. This correspond to the fact already 
noticed that it is only the global evolution of the interface and not of its parametrization that matters. 
Different choices of the parametrization lead to the following formulas: 

• When a = (this corresponds to p- = p+ i.e. to the Kelvin- Hclmholtz problem), it is convenient 
to choose a Lagrangian parametrization, 

r t (X, t) = v(r(X, t), t) = U+ ± U - (r(A, t))\ (27) 

the equations (|23|) and (|2^|) then reduce to 

d t u = and thus w(A, t) = w(A, 0). (28) 

Using p5)l . we thus get 

<V(A, t) = p.v. I V^G(r(A, i), r(A', t))«(A', 0)dA' (29) 

or, when Q = M 2 , 



ftr(A, t) = ±-R« p.v. / ,;A'; ; Z7^a "(^» °) dA ' • ( 30 ) 



1 f r(X,t)-r(X',t) _, i; 

2^ fP - V 7 |r(A,t)-r(A',t)| 

It is therefore natural to reparametrize the interface by the (time independent) arc length of its 
initial position, T,(t) — {r(s ~ 1 (s),£)}, where sq(X) = f Q d\r(X',Q)d\'. When 51 = R 2 , one can 
further identify the fluid domain with the complex plane, i.e. z(s,t) ~ t(sq {s),t) to obtain: 

d t z(s, t) = -Lp.v. / ll S '' 0) ds', (31) 

27T2 J 2(s,i) — Z(S , t) 

where j(s, 0) is the vortex strength (IT21 at t = 0. 

If moreover 7(s, 0) has a distinguished sign (which is the case when the initial vorticity is a Radon 

measure), then we can use the mapping 

ao0)= / 7(s',0)ds' 
Jo 



to perform another (time independent) reparametrization of the interface, namely, E(t) = {z(a (a),t)}. 
One thus obtains the Birkhoff-Rott equation for the variable z(a, t) = z(a^ (a),t), 

^"-^■/ zMl-^V) - (32) 

The Birkhoff-Rott equation corresponds to a parametrization by the circulation; indeed, one readily 
checks that \d a z(a a (s),t)\ = ^gy- 

• When a = — 1 (this corresponds to p + = 0, and thus to the water waves problem), several choices 
are possible. One can for instance, as in [54] . choose a Lagrangian parametrization with the velocity 
given by the trace of the velocity field in Sl—{r) at the surface, 

r t (X,t) =«_(r(A,t),t) =u™{r(\,t),i)n + v?L(r(\,t),t)T. (33) 

Another choice, adopted in [5], consists in parametrizing the surface by its (renormalized) arclength. 
It is also possible to use the particular case a = — 1 of the following point when the surface is a 
graph. 

• When the interface is a graph (we have seen that this hypothesis can always be used for "local 
analysis"), we choose the canonical parametrization by x\, S(t) = {(x\,<j{xi,t)} and the full system 
(including the Biot-Savart equation) becomes: 

d t a + vid Xl a = v 2 (34) 

d t (-uj - a(vi +v 2 d Xl a)) 

1 ~ 2 I |2 

+d Xl |wi (-2 - a{vi + v 2 d Xl cr)) - a ( 8(1 + g a 2\ ~ — ~ giXl ~ • 92Cr ) } = ° ' ^ 35 - ) 
/ .x 1 f cr(xi,t) - a(x[,t) , , 

v ^) = -^-j ixi ^ x[) 2 + {aixut )^a( x ' 1 ,t)y UJ ^ t)dx ^ (36) 

4 Linearization of the interface problem 

To analyze the local properties of the above system (l34l) - (|37p . one introduces its linearized version near a 
stationary solution a(x,t) — 0, u,(x, 0) = (— uo,0), and u_(x,0) — (uq,Q) (and thus w = 2uo)- Looking 
for solutions of the form a + eo~, w + eu5, we are led to the following equations: 

3 t o- - v 2 = , (38) 

d t (-ui- avi) + d xi (uqV! - -au u + ag 2 a) = , (39) 

evi(a:i,t) = -e^p.v. / a ^ Xl ' t) - ^f 1 '^ (2u Q + eui(x' 1 ,t))dx' 1 = -eu Q \D\a + Q(e 2 ) , (40) 
27T 7 (xi-x'J^ 

ev 2 (x 1 ,t) = — p.v. f J^—±- (2u + ew(ri,t))da:i = ^e-ffw + 0(e 2 ), (41) 

(42) 



where the operators \D\ and H are as defined in ([22)) and (|2ltf (and where the identity |D|1 = HI = 
has been used to derive the second part of (|40)l and ([H])). 

At leading order in e, we thus obtain the two following equations: 

d t a-~HZ = 0, (43) 

d t (-uj + au \D\a) +d Xl ( — u^\D\a - -cluqlj + ag 2 a) = 0. (44) 

Taking the derivative with respect to t in (|4"4l and using (|4U)) one obtains finally 

d 2 u - 2au Q d 2 Xi fi + u\d 2 x Z) + ag 2 \D\Zj = 

or 

{d t - au d Xl ) 2 uJ + u 2 {\ - a 2 )dl x u + ag 2 \D\u = . (45) 

The nature of Equation (|45[) depends on the parameters a, g 2 and uq: 

• If the Atwood number satisfies a 2 < 1 and if u 7^ 0, then the equation (J45I) is always elliptic, 
regardless of the sign of ag 2 ■ 

• If a 2 — 1 and ag 2 > 0, which corresponds to the water waves case (with downwards gravity), the 
equation (I45|) is then nonstrictly hyperbolic and the condition ag 2 > is in fact the Levy condition 
ensuring the wellposedness of the Cauchy problem associated to it. 

• If «o = and ag 2 > 0, which corresponds to the Rayleigh- Taylor case (0 < a < 1) with the heavier 
fluid placed below the lighter one, then (j45|) has the same characteristics as in the case a 2 = 1, 
ag 2 > 0. The ill-posedness of the Rayleigh- Taylor problem, even for configuration close to the rest 
state, are therefore due to nonlinear effects. 

• If a 2 = 1 (or «o = 0) and ag 2 < 0, which corresponds to the water waves case (or Rayleigh- Taylor) 
with the lighter fluid below the heavier one, the equation (|45[) is then hypoelliptic. 

Remark 2. A Fourier type analysis of (|45[) with the introduction of modes 

uj(x,t) = e Tt+i ^ 

shows that all frequencies £ such that 

Igl > Qg2 

lead to unstable modes (with 5Rt > 0) . The range of unstable frequencies increases with a going to zero 
which is the Kelvin-Helmholtz case. Without the nonlinear interactions, low frequencies would therefore 
remain stable for the Rayleigh- Taylor problem. 

5 Consequences of the ellipticity of the Rayleigh- Taylor and 
Kelvin-Helmholtz problems 

5.1 Analytic solutions 

For a precise statement of the behavior of analytic solutions to (l34l) - ([37|) . it is convenient to follow 
Duchon and Robert |19j and to introduce B p , the space of analytic functions which are Fourier transform 

10 



of bounded measures with exponential decay at infinity (keeping in mind the Payley Wiener theorem, 
this corresponds to functions that can be analytically extended in the strip {(x + iy) , \y\ < p}) : 

B P = {u = f d\u\, \\u\\ p = f d\e p ^u(£_)\ <oo}. (46) 

One then has the following theorem for the Kelvin-Hclmholtz and Rayleigh- Taylor problems; it shows 
that one can construct local in time solutions to (|3~4"1) - ([3"T)) if both the interface and the vorticity are 
prescribed (and analytic) at t = 0, and that global in time solutions exist if the initial vorticity is not 
arbitrarily chosen. Note that in the second point of the theorem, one considers a system on (d Xl a, w) by 
taking the x\ derivative of ((34)) . 



Theorem 5. 1- (see \5S§ and \5J$ ) Let p > and (cot^o) £ B p , and assume further that uq has zero 
mean value (i.e. (To({0}) = 0). Then there exists a constant a such that the Cauchy problem J^p - f^T} ) 
has a unique solution 

(a(t,.),ui(t,.)) eB ap - W 

defined for \t\ < ap (and whose dependence on t is analytic). 

2- (see U9f ) There exists a constant e > such that for all d xi o~o € Bq satisfying 

/"|» (OI<e and C^({0}) = 0, (47) 

there exists a > such that the problem d Xl J34\ )- Ij37\ ) admits a unique solution (d Xl a,uJ) € C(0,oo;£> Q ) 2 
with the boundary conditions: 

d Xl a(0,xi)=d xl a (xi), S5(0p)({0}) = 0. 

Moreover, d Xl o~{-,t) and w(-,t) converge to uniformly as t goes to oo. 

Remark 3. 

1- The first point has been proven first for the Kelvin-Helmholtz problem (a = 0) in [55] and then extended 
to a ^ in [541 . The basic new ingredient is contained in the relation between Co and u/2 — a(vi+V2d Xl o~) . 
The term (v± + V2d Xl a) is proportional to v T the tangent component of v and one has: 

(vi +v 2 d Xl a) = K(lo) 

with K being a regularizing operator as observed in Remark [TJ Therefore with Fredholm property and 
spectral analysis one can, for \a\ < 1, invert the relation 

-ui - a(v! + v 2 d Xl o-) = (-I-aK)(ui) (48) 

to obtain a well posed Cauchy problem for (cr,ui) in the class of analytic functions (see for details [28] . 
[54] and [8]). 

2- The proof of the second point has only been done for a = but the extension to < \a\ < 1 does not 
seem too difficult. 



If 



Looking at the solution provided by the second point of the theorem on the interval [0,T], with 
T > 0, and changing the time variable t *— > T — t, we have at our disposal solutions which are analytic 
for < t < T and which at the time T exhibit any type of singularity compatible with the hypothesis 
(|47p . This is a much larger class of singularities than the one initially predicted by Meiron, Baker and 
Orsag for Kelvin- Helmholtz [44] . One can have examples for which the Hausdorff dimension of the graph 
%i ^ d Xl a{x\, T) may be of dimension greater than one. In particular in a periodic version of the above 
theorem one obtains starting, from analytic initial data at time t = 0, a solution that coincides at t = T 
with the Weierstrass function, 

a(x u T)= J2 2-^~" cos^zx), ft^&T) = £ 2"* (*(£ + 2") + 8(t - 2")). 

0<n<oo 0<n<QO 

The computations made in [U] and [33] led to the belief that the first appearance of singularity was 
the formation of a cusp in the curve d Xl a(xi,T) or u>(x\,T). As already said, this turns out to be a 
special case of the above theorem, but that was also the motivation of the construction of Caflisch and 
Orellana [TU] ■ With perturbation methods very close to the one used in [TS], they construct solutions 
for Kelvin-Helmholtz problem which are analytic for t > and which at the time t = have an isolated 
singularity behaving like a cusp. In complex representation these solutions are constructed as analytic 
perturbation of the function: 

z (A, t) = (1 - t){(l - e-i- jA ) 1+ " - (1 - e-3+ iA ) 1 + 1 '}, (49) 

which turns out to be an exact solution of the linearized Birkhoff-Rott equation and which has the 
following properties: i) For any t>0Ai-} Zo(t, A) is analytic, ii) For t = A h- > z(0, A) does not belong 
to the space C 1+u but belongs to any space C 1+v with < v' < v . 

5.2 After the appearance of the first singularity 

The mathematical properties of the interfaces may be used (with a contradiction argument) to obtain 
some results concerning the behavior of the solution (if it exists) after the appearance of the first singu- 
larities. As this turns out to be in many free boundaries problems, the method consists in showing that if 
the solution is locally (in space time) more regular than some threshold 5 then it is very regular (locally 
C°° or analytic). This implies in turn that after the appearance of the first singularity, the solution 
has to be less regular than the threshold S. As a consequence, the determination of a threshold S that 
would imply C°° or analyticity is one of the challenges of the theory, and as such has been studied for 
instance by S. Wu [60] , Lebeau [33] and Lebeau-Kamotski [28]. For instance in [28] one has the following 
statement: 

Theorem 6. Assume that locally near a point (so,io) the curve S(4) = {r(s,t)\ parametrized by the 
arc length belongs to C 1+l/ (v > 0) and is chord- arc; assume also that the density of vorticity is strictly 
positive and bounded at (so,£o)- Then in the neighborhood o/(s ,io); the curve and the density of vorticity 
are C°° . 

Remark 4. If the vorticity is assumed to be positive and bounded for all s (and not only in a neighborhood 
of so), then the curve and density of vorticity are analytic. It is expected that the conclusion of Theorem 
[6] remains true if one replaces C°° -regularity by analyticity. 
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The result of S. Wu [SD] obtained independently and previously for Kelvin- Helmholtz is slightly 
sharper. Instead of C 1+l/ one assumes that the function r(t, s) belongs to H 1 and does not rolls up too 
fast. More precisely one has the following theorem (we give a non sharp version of the condition (J5H)) to 
avoid the definition of technical functional spaces). 

Theorem 7. Assume that z £ ff ^[0, T], Ifo c (R)) n L 2 ([0, T]; ff/ M (l)) is a solution of the Birkhoff-Rott 
equation liSS\) , which is uniformly chord-arc on [0,T]: 

3m,M>0, V(i,a,/3) e [0,T] x M 2 , m\a - j3\ < |z(a,t) - z(/3,i)| < M\a - /3|. 

TTien i/iere exists a constant c(m, M) such that the relation: 

sup | ln(z Q (-, t)) |U <c{m,M) (50) 

[0,T] 

implies the analyticity. 

Remark 5. The above statements are in full agreement with previous theoretical and numerical results. 
The example of Caffish and Orellana corresponds to a singularity in C 1+lJ and the Lebeau theorem shows 
that this function cannot be extended as a solution in C 1+v with v' > after the singularity. 
The Prandlt-Munk vortex, described in Section [5J below, is a non regular solution of the Birkhoff-Rott 
equation defined for all time. This is not in contradiction with the previous theorems because near the 
points (±1, |) the vorticity density does not satisfy the Lebeau- Wu hypothesis. Numerical simulations 
done by Krasny [30] indicate the appearance after the singularity of spirals. On the other hand Theorem 
[7J implies that some spiral- type behavior would imply regularity and therefore is excluded. The catch lies 
in the fact that Krasny spirals do not seem to satisfy an hypothesis of the type (|50l) 

Without entering into details it is interesting to give some comments on the proof of Theorem UJ 
As in any free boundary problem the challenge is to determine the threshold of regularity S as sharply 
as possible. The most delicate step is to prove that this threshold regularity implies a slightly better 
regularity; deducing analiticity (or C°° regularity) from this slightly better regularity is usually much 
easier. 

In the present case, one shows that the assumptions of Theorem [Jj lead to iJ 2 -regularity. This can 
be reached from the Birkhoff-Rott equation by a localized (with a smooth real valued function rj(a)) 
variational estimate. With V — i](a) (1 + i) lnz a one obtains the formulas 

V t-^^\D\V^F(a,t), 

Z|Zq,| 

for some function F(a,t) that can be explicitly computed, and thus 

/ |Zq|2| ( F * " ^p \D\V) 2 da = J \z a \ 2 \F(a 7 t)\ 2 da. (51) 

On the other hand, one easily checks that 

J\^\ 2 \( V t-^J 2 \DW) 2 da = J (\ Za f\ Vt \ 2 + ^\\D\V\ 2 ) 
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Since the assumptions of Theorem [7] imply that m < \z a \ < M, it is possible to deduce (after a bit of 
work!) from (15T1) and (|52|) , a control of the quantity 

q (t) J ( m z\V t f + ^\\D\V\ 2 )dadt, 

where q(t) is a carefully chosen smooth, real valued function. Recalling that V = r](a)(l + i) lnz Q , this 
suggests that it is possible to control the -ff^-norm (in space and time) of lnz Q by its L°°-norm, which 
is finite by assumption. 

6 Nonregular vortex sheets 

6.1 The Prandtl-Munk vortex sheet 

As described in Section [3j a solution of the Kelvin-Helmholtz problem (computed on the interface) 
should provide, with u given by the Biot-Savart relation (|14j) . a weak solution of the incompressible 
Euler equation. The situation turns out not to be so simple for the following reasons: 

(i) Relaxing regularity hypothesis on the interface, one can generate solutions of (|3T)|) which are not 
solutions of the Euler equation. And this pathology appears with initial data which satisfy the 
hypothesis of [39] and therefore have weak solutions (different from what is given by the Kelvin- 
Helmholtz equation ([50])'). 

(ii) Even without the above pathology, one should observe that the Kelvin-Helmholtz problem generates 
solutions which do not satisfy (and there is a big "gap" ) the uniqueness stability criteria of Yudovich. 

The point (i) is illustrated by the "Prandtl-Munk vortex sheet" which is notorious in fluid dynamics as 
the one that generates a circulation distribution that minimizes the induced drag on a plane wing |46) . 
With initial data: 

E(0) = Mn.O) = (si,0)}, w(a:i,0) = 7 JL= X (-i,i)(si), 
where X(-i,i) is the characteristic function of the interval (— 1, 1). One observes (cf. [46]) that 

r(xi,t) = (xi,--), v = (0,--), w(xi,i) = u(xi,0) 

is a solution of the equation (|30|) . However the corresponding velocity u(x,y,t) is not a solution of the 
Euler equation (fTUf and one has instead (see [55]): 

V • v = 0, d t u + V • (u (g) u) + Vp = F, 

with a nonzero forcing term F given by (with 5 denoting the Dirac mass): 

F(xi,x 2 ) = - ((*(xi+l,» a +4) ^ *(»i-i,* a +4)). °) ■ 

The catch lies in the fact that the vorticity density is not regular enough in the neighborhood of 
the points x\ = ±1 to justify the formal computations. On the other hand the initial data satisfy the 
hypothesis of [391 and there is therefore a weak solution of the Euler equations with these initial condition. 
This solution is singular enough near the same points to differ from the one given by \30\) . We refer 
to [37] for a discussion of the equivalence at low regularity of the Birkhoff-Rott and Euler description of 
vortex sheets. 
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6.2 Regularization of vortex sheets 

There exist various ways to regularize a vortex sheet, with different properties. 

(i) Convolution by a sequence of functions converging to the Dirac distribution of analytic initial data 
with vorticity concentrated on a curve provides a regular solution of the Euler equation that does 
converge to the corresponding solution of the Kelvin-Hclmholtz problem (in the existence domain 
of such solution) , see [JT] . 

(ii) The so-called a-regularization (one of the novel approaches for subgrid scale modeling of turbu- 
lence): In R 2 the vorticity transport equation: is replaced by the a regularised equation: 

d t uj Q + u Q ■ \7uj q = Q,u a = K a u) a = (I-aA) _1 M with u(x,t) = — -% / , pruj a (t,y)dy . (53) 

2tt 2 J R 2 \x - vY 

For any a > the kernel K a is regularizing. Therefore the above system is well posed (with global 
in time existence and uniqueness of the solution) for any measure valued initial data. Moreover 
any initial vorticity density supported by a curve S(0) remains on the curve T, a (t) transported by 
the flow 

X(x, t) = u a (X(x, t),t), X(x, t = 0) = x. 

Eventually if the initial data (curve and vorticity density) are analytic, for a — >• the curve S a (i) 
and the corresponding vorticity density converge to the solution of the Kelvin Helmholtz problem 
as long as this solution, which exists locally in time, is defined; for details and further references 
see 0. 

(iii) However and more annoying is the following fact: with an initial vorticity concentrated on a curve 
the Navier-Stokes has with viscosity v > a well defined solution u v (x,t). Furthermore if the 
initial data satisfy the hypothesis of [39} this solution converges to a weak solution. 
Weather this solution coincides with the one given by the Kelvin-Hclmholtz problem (even with 
convenient hypothesis on the regularity of the curve and of the density of vorticity) is a completely 
open issue. With the "Prandtl-Munk vortex sheet" one can exhibit a counter-example (but with 
singular curve and vorticity) to this convergence: The corresponding initial data satisfy the hypoth- 
esis of [35] • Therefore, when the viscosity goes to 0, the corresponding solution of the Navier-Stokes 
converges to a weak solution of the Euler equation. On the other hand the "Prandtl-Munk vortex 
sheet" which is the solution of the Kelvin Helmholtz equation with the same data is not a weak 
solution of the Euler equation. 

7 Water waves 

The water waves problems refers to the case (like water in air) where one of the density (say p+) 
is infinitely small with respect to the other and then taken equal to zero. As said in A3. 21 various 
parametrizations are possible. Parametrizing the surface £ = {r(t, A)} as in (|33|) . one has 

a = — 1, r t — u'ln + u t _t, r\ = f\r\\ . 

Therefore in the system ([2"5]) - (|2^)) . or equivalently 

(d t r — v) ■ ft = 0, 

15 



one may use the following relations 

Cj 

- + {v-rx) =u r _\rx\, 

1 / n ,& , .. ,1 _ _ \u T _\ 2 . 
r-p (u - r 4 ) ■ r A (- + (v • r x )) = {-u_u + — ) , 

l8|r A | 2 2T l 2 M -"+ 2 j ' 

a t K|a A r|)-9 A (i(| U n | 2 + |^| 2 +g.r)=0. (54) 

A variant of (f54| can also be derived if one chooses to parametrize the surface by its arclength (see for 
instance [22]); it is also worth commenting on the case where the surface is a graph, £(£) = {{x\, c(x\, t))}. 
Starting from (|34p -([35 |) and remarking that 



to obtain the equation 



-w + (vt + v 2 d xi a) = \/l + {d xl a) 2 u T _, 
one gets 

dt(y/l + (d Xl a) 2 uZ_) + d Xl [u r _(v - r t ) ■ f } + d xi {±(u r + - u T _) 2 - ±\v\ 2 } - g ■ r Xl = 0. 
Since rt = (0, dtu) T — (0, \/l + (c^er) 2 ?!™ ) T , this is equivalent to 

a t (\/l + (& lf r) 2 U l)+a a;i {i( U :) 2 - t£t£ - i(,/!) 2 +. g( r} = 0, (55) 

where we assumed that g = (0, —g) T - 

Since V • u = and V A u — in fi_(i), and assuming that Q_(t) is simply connected, there exists a 

velocity potential 3>(x,£) which satisfies the relation: 

V$(z,i) = u(x,t), A®(x,t)=Q in fi_(t). 

Denoting by 'J the trace of <£ at the surface, \fr(a:i,£) = &(xi,<r(xi,t),t), one has 



On* = a xi $| E(t) + 9 Xl cr9 X2 $| E(t) = v/l + ^a) 2 ^. 
Integrating with respect to Xi in f|55[) yields therefore 



dt$ + l( u -) 2 ~ u -u 7 l - \( u -f + 9° = ° (56) 

(note that adding if necessary a function depending on t only to $, the integration constant can be taken 
equal to 0). 

Eventually it is convenient to introduce the Dirichlet Neumann operator G(a) : H's(M.) n- H~z(M.) 
associated to the resolution of the Dirichlet problem in f2_(£) C f2 (with an additional homogeneous 
Neumann condition corresponding to (j4]) if there is a fixed boundary) : 



G(<r)tf = y/\ + (d Xl a) 2 d H $, (57) 
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where 



Remarking that 



-A$ = in Q(t), 

$ = * on dQ-(t)\{dn_(t) DdQ), dn® = on (0fi_(t) ndil). 



y/l + (d xx <r) 2 u T _ = d xl V and y/l + (d mi a)H» = G(a)*, 
we deduce from (J56I) that the water waves equations can we written 

W + ga+ l -{d x ^f - 2(i + (0^)2) ( G ( CT )* + ^^xx*) 2 = 0, ' "'' ; ! 

where we assumed that g = (0, —g) in the above expression. 

Remark 6. The formulation (I58[) of the water waves equation corresponds to the Hamiltonian exhibited 
by Zakharov |63j and was written first under this form by Craig and Sulem [16] . As it is well known the 
equation (|58[) can be obtained directly from the Euler equation 

dtu + u-Vu-i Vp = g (59) 

written on the domain f2_(£) occupied by the fluid. The above short derivation has been written to com- 
pare this problem with the Ray leigh- Taylor and Kelvin- Hclmholtz equations and to emphasize similarities 
and differences. 

7.1 Rayleigh- Taylor type behavior 

In the water waves case, the linearized equation (1451) . with uq — 0, can be written 

d^u~g 2 \D\u = Q. (60) 

It has already been noted in JjHthat if 52 > (the fluid is above the surface), the operator appearing 
in (ISUl) is hypoelliptic. The observations of Section [5] involving the ellipticity of the Kelvin Hclmholtz 
problem could be adapted. To the best of our knowledge this has not yet been worked out with full details 
except in the case of traveling waves which is the object of §7.^[ and ^8.2\ for the extension to the 3d case. 
Near a stationary solution with g 2 > one may prove, with analytic initial data, local in time existence 
of an analytic solution. Relaxing one of the constraints at t = one may adapt the construction of |19] 
or [TU]. Hence by the change of variable t M> T — t, one could construct solutions that exhibit a large class 
of singular behaviors. Finally since a Galilean transform changes the problem into an equation involving 
perturbations of a linear term of the form (|60|). the results of [33] and [60] could be adapted leading 
to a threshold of regularity S which for solutions more regular than this threshold would have C°° or 
"Gevrey" regularity. In some sense the water waves problem with agi < shares the same behavior as 
the Rayleigh- Taylor problem and therefore carries the same name. 

7.2 The nonstrictly hyperbolic behavior 

When gi < (the fluid is below the surface), we have already commented in g]that (ISUl) is nonstrictly 
hyperbolic, and that the condition gi < corresponds to the Levy condition on the subprincipal symbol. 
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The equation (|60|) roughly behaves as a wave equation with the positive self-adjoint operator —&%. replaced 
by its square root. 

An important observation is that the sign of gi is an intrinsic object. The equation (1601) has been 
obtained by small perturbation analysis near a point where the x\ axis is tangent to the interface and 
where the interface behaves like (xi,a(xi,t)) therefore gi coincides (modulo small perturbations) with 
g-n. With the full, nonlinear equations (|5§)) . the corresponding quantity is —dnp. Therefore the condition 

- dnp > on S(t) (61) 

should enforce the stability of the problem. This condition is called the Rayleigh- Taylor stability criterion 
(or Taylor sign condition). From a mathematical viewpoint, it is also the Levy condition ensuring the 
wellposedness of a nonstrictly hyperbolic equation. This is transparent on the formulation of the water 
waves problem given in Proposition [5] below. Following [21], we denote by D t = dt + v ■ V the lagrangian 
derivative, and by M the operator A/" = Wl + (d Xl a) 2 G(cr) (note that it is easy to define M even if £(£) 
is not a graph); one then proves the following result for the interface curvature k: 



Proposition 8. (cf Proposition 2.6 in \2J^] ). For any smooth interface solution of the water wave 
equation one has: 

D 2 t n + (-d H p)Afn = R, (62) 

where R is a remainder term consisting of lower order terms. 
To prove this proposition one considers 

d t u_ + u_ ■ Vu_ H Vp- = g (63) 

9- 

on the interface dfl^(t) and apply the Laplace Beltrami operator Ajv t ) to (f6"3"|) ; using the important 
remark that at leading order, the Dirichlet-Neumann operator G{a) coincides with the square root of 
the Laplace Beltrami, 

G(a) = J-As(() + Ro, Ro — zero-th order operator, 

one can establish that modulo lower order terms one has: 

ft- A E ( t )(Vp(t,a:(t,A)) =dnpNn, 

and the result follows quite easily. 

With the formula (|62l) one obtains local in time stability of the interface for initial data having a 
limited Sobolev regularity. Of course, this result can be obtained with other techniques (see for instance 

[sg GU ED nana]). 

Remark 7. In the whole space R the solution of the model dispersive equation (corresponding to (|60[) 
with (72 = — 1) 

dfu+ \D\u = } u(x,0)=0, u t (x, 0) = <j>(x), 

satisfies a dispersive estimate 

\\u(;t)U~<Ct-i\\<j>(x)\\ L r. (64) 
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Therefore, one may expect to use dispersive effects, when the interface is a graph, to obtain large time 
estimates for the solution. This has been done in 2d (i.e. surface dimension d— 1 = 1) by S. Wu [61 . The 
central idea in [5T] is to transform the equations (inspired by Birkhoff normal forms) in order to increase 
the degree in the non linear terms and therefore to reduce their size. Though considerable difficulties 
make the implementation of this programme very delicate, it leads to an existence time of the order of 



exp 



C 



some norm of the initial data 



Since the dispersive effect in the formula (|64[) improves with the dimension, the situation for this problem 
is better in 3d (i.e. surface dimension d — 1 = 2, see 



7.3 The Ray leigh- Taylor stability criterion 

Our goal here is to comment on whether the Rayleigh- Taylor stability criterion (|61|) is satisfied for the 
water waves problem. Let us consider the case of water waves in finite depth. The boundary of Q- (t) is 
the union of two distinct pieces. A given boundary Tfi X where the fluid is in contact with the boundary 
of the vessel, and the free boundary E(i) between the region of the fluid of density pi (we take p- — 1 for 
the sake of simplicity), and the fluid of density p + = 0. On Tfi x one assumes as usual the impermeability 
condition. One thus gets 



d t u + V(u u) = -Vp + g, V-u = 0, VAu = in Q-(t), 

p = on S(i), u-n = on Tfix- 



(65) 
(66) 



Note that neither £(£) nor V fix are assumed to be "globally " a graph. However one consider only 
configurations where 

d(T fix , E(«)) > h m in > 0, (67) 

and therefore where the influence of the fixed boundary on the free boundary involves only lower order 
terms (in terms of regularity) . 
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Taking in particular into account the relation V A u_ = 0, one deduces from (|6"5)) - (l6"6l that the 
pressure p satisfies 



-Ap = J2a d Xi Ujd Xj Ui = |Vw| 2 > in fi_(i), 



p = on S(i), dftp = — /lujUj-z-^ + g • ^ on 



fix- 



*3 

Then the maximum principle implies that p > in 0_(t) and therefore that 

-<9 s p>0 on £(i), 
as soon as 

f)~n = \ 11:11 

dx 



«sP = ~ > J u i u jj^r + S-n > on T fix . (68) 



For instance (l68l) turns out to be true for flat bottoms, since one then has 



Edfii „ ^ _ ^ 
UiM i "5 hg-n = g-n>0; 
ax, 
y 

the same argument works also in infinite depth (in fact, this is the configuration where it has been proved 
first [57J[SH]), or for liquid dropH 

Remark 8. Flat, infinite bottoms, or liquid drops are of course not the only configurations where the 
Rayleigh- Taylor stability criterion ([BTjl is satisfied. Quite obviously, the same arguments show that (fBTj) 
is also satisfied for "almost" flat bottoms. More generally, a condition on the curvature of the bottom 
can be derived [23 [3] ■ 

7.4 Ellipticity and regularity of traveling waves 

Ellipticity results have been used above to derive instability results for the Kelvin-Helmholtz and Rayleigh- 
Taylor problems. Though such instabilities do not occur for surface waves (when the Taylor sign condition 
(|6ip is satisfied), an elliptic behavior can still be observed in some phenomena concerning water-waves. 
In particular, it is possible to derive a priori regularity results for traveling waves, both in 2d and in 3d. 
Comments on the 3d case are postponed to ij8.21 

Using the formulation (1551) of the water-waves in terms of the surface elevation a and the trace of the 
velocity potential at the free surface ^, a traveling wave is a solution to ([58"]) of the form 

(a(t,xi),ip(t,Xi)) = (a(xi - ct),-ip(xi - ct)), 

where cGlis the speed of the traveling wave. Consequently, (a, ip) must solve (we omit the underlines 
for the sake of clarity) , 

G(a)ip + cd xi a = 0, 



3 However, this argument fails if the fluid is not irrotational. This is the reason why, in the study of liquid drops, it is 
assumed in 36 , 49 , 13 that the Rayleigh- Taylor condition holds. It also known that without this condition, the equations 
are ill-posed 1201 
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Existence of traveling waves in one surface dimension is known from the early works of Levi-Civita 

[34] and Nekrasov [45] . The fact that if the free boundary of these solutions is a C 1 curve, then it is also 

analytic is also a classical result that goes back to the works of Lewy [35] and Gerber [2U [33J [23] . These 

results use the fact that the system (IM|) is elliptic, a fact that can easily be observed on the linearized 

equations around the rest state a = 0, ip = 0. Recalling that G[0] — \D\, the linearized equations read 

indeed 

\D\ip + cd Xl <7 = 0, 

go - cd Xl ip = 0. 



The symbol of this system is therefore 



9 -ic£i 



whose determinant is the elliptic symbol c £% — g\£ 
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7.5 A word on shallow water asymptotics 



Except for its influence, through ([68]). on the Ray leigh- Taylor criterion ([6l]) . the bottom does not play 
any role in the well-posedness theory of the water waves problem described above. By assuming (|67[) . i.e. 
that the water depth never vanishes, one is ensured that the contribution of the bottom to the Dirichlet- 
Neumann operator G(cr) 1 I' is analytic (by the ellipticity of the potential equation A$ = 0); therefore, the 
contribution of the bottom is ignored in any symbolic description of G(a)^ . This is the reason why it is 
possible to obtain a local well-posedness theory for very rough bottoms, as shown recently in Q]. 

For the analysis of instabilities that interests us here, the influence of the bottom is thus irrelevant. 
This fact seems to contradict common sense observations, which show that the bottom often plays a 
important role in the dynamics of the surface (consider the shoaling of a wave as it approaches the shore 
for instance). The explanation of this apparent paradox is that the fact that the contribution of the 
bottom is smooth (and can therefore be ignored in the wellposedness theory) does not mean that it is 
small. For instance, consider the approximation of the Dirichlet-Neumann operator by the square root 
of the Laplace Beltrami of the surface, 



G(a) = </-Aj;([) + Ro, Rq = zero-th order operator; 



we have already seen that this symbolic approximation is very useful to prove the stability of the water 
waves problem. However, it is not uniform in the shallow water limit. More precisely, the operator 
norm ||-RoIIl 2 ^l 2 grows to infinity in the shallow water limit (which corresponds to configurations where 
the ratio of the depth over the typical wavelength of the interface deformation is small). In order to 
obtain uniform estimates in the shallow water limit, one must take into account the contribution of the 
topograph^]. In some sense, the study of the interface is not decorrelated from the topography in the 
shallow water regime (the same phenomenon also occurs for interface waves, see for instance [llj). In 
order to focus on the instability phenomena at the interface, we implicitly assumed throughout this paper 
that we are not in a shallow water type regime. The study of this latter is however possible, albeit 
with different techniques, and a consequent literature has addressed this problem. A few examples are 

na i2si eh 123 ei . 



4 In 1321 , a symbolic analysis "with tail" of the Dirichlet-Neumann operator is proposed, the tail being an infinitely 
smoothing perturbation of the principal symbol ./— Aw t \ that takes into account the contribution of the bottom and 
removes the shallow water singularity 
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8 Differences between the 2d the 3d case 



Most of this review has been devoted to interfaces in planar flows. Before concluding and without going 
to the details one should single out differences and similarities between the 2d and the 3d case. 
In 3<i the formula for the propagation of the vorticity is 



d t w 



\7uj 



Vw; 



as a consequence the 2d conservation relation does not exist any more. Therefore, there is no theorem for 
weak solutions nor global in time result for strong solutions. The only general available result concerns 
local existence in time of smooth solution, say, in C 1,a . This regularity seems optimal as shown by the 
instability of the shear flow [40 [ [3] : more precisely, consider the vector field 



U sh f(x,t) = (ui(x 2 ),0,u 3 (xi -tui(x 2 )), 



(70) 



which is an obvious solution of the 3d incompressible Euler equation. With this type of flow one can 
construct solutions with the following property 



t/ shf (x,0)eC 



0,o 



and for any t > 0, 



L/smOM) ^ C°< CT . 



8.1 About the 3d Kelvin-Helmholtz and Rayleigh- Taylor problems 



For smooth solutions, the generalization of the formulas (|25p - ([26|) for the equations of the interface remain 
valid. Hence one can prove a local in time existence theorem for analytic initial data (cf. [54]). 

However the ellipticity properties cease to be valid. As previously observed, these properties follow 
from the fact that the equation can be locally considered as a perturbation of the linearized problem. 
This linearized problem is elliptic in 2d. In 3d, one can still reduce the problem (see [55J HI]) to the 
situation where the interface is a graph 23 = a(x%, X 2 ,t), and where the solutions are perturbations of 
the stationary state a = 0, QJ.(xi,x 2 ) — (£.i,w 2 ,0). Then in Fourier variables the linearized problem in 
the Kelvin-Helmholtz case (the Rayleigh- Taylor case can be studied in the same way) is 



\ 



/ a 

UJ 2 
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where the operator A(£) is given by (denoting £ = |£|(cos0, sin0)) 
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■ sine 



-!l£| 2 |w| 2 sin0 
il£| 2 |w| 2 cost? 
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\ ' -§(£.w)sin( 

The eigenvalues of the matrix A(£) are 

(with multiplicity 2), 



I (£ • S) cos t 



i(£-£J)sin0 

f (£ "£) cost9 



(71) 



and 



±2^ A ^I 



Therefore, the operator dt — A(£) is no more elliptic and perturbations transverse to the direction of 
streaming are thus unaffected [12] , 



22 



8.2 About 3d water waves 

8.2.1 A global existence result in 3d 

Local existence for 3d (surface dimension d— 1 = 2) water waves under the Rayleigh- Taylor criterion (|5Tj) 
can be proved with tools that do not depend on the dimension [551 [33 HS1 EH] ■ This is of course not the 
case if one wants to prove a global existence result based on the dispersive properties of the equation. 
In order to describe the stability result for the water waves in 2d (surface dimension d= 1), we used the 
equation on the curvature k provided by Proposition [8] In [24], it is proven that such an equation still 
holds in 3d with the curvature replaced by the mean curvature, and a careful analysis of this equation 
allows the authors to obtain a global existence result for small datas. This is achieved by using several 
ingredients (in the same spirit as those in 2d, see Remark [3 but with different techniques): 

• In the whole space R 2 , the solution of the model dispersive equation 

<9 2 u + \D\u = 0,u{x,Q) =0,u t (a;,0) = 4>{x) , 



satisfies a decay estimate 



\\u{x.,t)\\L~<Ct- l U\\ L i. (72) 



• 



The use of a Birkhoff normal form to increase the degree of the non linearity. However, time 
resonances do not allow the construction of a standard normal form; the authors therefore introduce 
the notion of space resonance, which catches the fact that, for dispersive equations, wave packets 
do not have the same group velocity and therefore "do not see each other" in general. 

8.2.2 Traveling waves in 3d 

Traveling waves also exist for 3d water waves (surface dimension d — 1 — 2), but this fact is considerably 
more difficult to establish than in the 2d case and has only been solved very recently by Iooss and 
Plotnikov |27j . As in the 2d case, a simple look at the linearized equations is very instructive. The 
symbol of the linearized system is indeed given by 

*c& 141 
9 -ic£,\ 

with |£| = (|£i | 2 + |^2 1 2 ) 1 / 2 - The determinant is then the symbol c 2 £ 2 — g\£\. This symbol is not elliptic 
since it vanishes on the (unbounded) region C — {c 2 £| = g|£|}. The equation defining C is a Schrodinger 
equation on the boundary that may propagate singularities. 

Despite of this, Alazard and Metivier [2] managed to prove that smooth enough biperiodic (i.e. 
periodic in both space directions) traveling waves are automatically C°°. By a careful paralinearization 
of the equations and two transformations inspired by the work of Iooss and Plotnikov, they show that 
the general situation is quite similar to what is observed in the linearization around the rest state. Then, 
they exhibit a diophantine condition that allow them to handle the degeneracy of the principal symbol. 

Remark 9. As already mentioned, the sign of gravity also plays an important role. Were it negative, the 
symbol of the linearized equation would be weakly elliptic, even in the 3d case. This remains true for 
the full equations (|69|) if the Rayleigh- Taylor condition is not satisfied. The paralinearization strategy of 
Alazard-Metivier can then be used. 
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9 Conclusion 

As a conclusion, we recall some of the basic threads that lead to the organization of this survey and add 
extra remarks and open problems. 

All the analysis presented here is deduced from the incompressible Euler equation written in the sense 
of distributions. With this setting the stability of the Cauchy problem has been analyzed in the models 
problems of Rayleigh- Taylor (Atwood number < \a\ < 1), Kelvin-Hclmholtz (a = 0), and water waves 

(H = i). 

With an high frequency analysis it has been shown that for \a\ < 1 the linearized problem is elliptic. 
Therefore the Cauchy problem is well posed only with analytic initial data. Moreover the solution 
whenever it exists is analytic even if the data (say for t = or T = t ) exhibit a singular behavior. To 
persist after the singularity the interface has to be very wild and this may be in agreement both with 
the numerical simulation of [30] and with physical observations. 

The situation for the water waves may be different. When the Rayleigh- Taylor sign condition —dftp > 
is not satisfied the linearized problem shares many similarities with the Rayleigh- Taylor equations. On 
the other hand with —dnp > on the interface stability results can be obtained. 

The above results given in 2d are compared in Section [8] with the situation in 3d. The ellipticity 
(or hypoellipticity) of the underlying linearized problem is lost and therefore the situation may be very 
different for the Kelvin-Helmholtz and Rayleigh- Taylor problems. As a companion problem observe that 
with the shear flow 

u(xt,X 2 ,X3,t) = (ui(x 2 ),0,u 3 (xi - tui(x 2 )), 

one can construct solutions defined for all time and with any "limited" degree of regularity on the 
interface. The same kind of difficulty is also encountered in some situation for the water waves problem 
(for the regularity of traveling waves for instance). 

In the present description physical quantities like the surface tension and the viscosity have been 
ignored. This is based on the observation that such quantities are smalQ and that the understanding of 
the pathology of the phenomena when they are made equal to should contribute to the understanding 
of what happens when they are small. 

• The introduction of a surface tension creates a third order dispersive term and all the above 
equation become locally in time well posed [26l HI [6] . This property does not persists when this 



5 this assertion is not very honest! Viscosity and surface tension are negligible for physical phenomena involving "rea- 
sonable" frequencies, but most of the instabilities described here occur at very high frequencies; the viscosity and surface 
tension terms should then play an important role. We refer to [32] , where it is shown that a natural two-fluid generalization 
of the Rayleigh- Taylor stability criterion 1610 is given by 



1 (P+P-) 2 
4T(p++p-)' 



-0*p]>t JV, '-^ WMi 



where T is the surface tension coefficient, u> is the jump of the horizontal velocity at the interface, and c(ct) is a geometrical 
constant. 

However this "not very honest mathematical behavior" is not restricted to fluid mechanics and nevertheless contributes 
to the understanding of nature. An other basic example is the Maxwell equations; they involve the permittivity e(x) and 
the magnetic permeability fi(x). These equations are the prototype of hyperbolic systems and the high frequency analysis 
is systematically used for their solutions. On the other hand e and fj, are derived from the Helmholtz equation under the 
hypothesis that the frequency is "moderate" . For higher frequencies k in the Helmholtz equation they would depend on k. 
Therefore in the space domain one obtains an equation with e and fi acting as time convolution and no more as coefficients 
and the nature of the equations at high frequency should be changed! 
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surface tension goes to zero except for the water waves problem with the Ray leigh- Taylor sign 
condition l5l l49l l42l . 



• With the presence of a viscosity v > the Euler equation is changed into the Navier-Stokcs 
equation. When the density is constant (p+ = p_), this equation is in 2d reasonably understood 
even for initial data with a vorticity being a measure concentrated on a curve. When v goes to 
zero (under reasonable hypothesis on the sign of the vorticity V A u v ) the corresponding solution 
converges to a weak solution of the Euler equation. Since there is no uniqueness theorem for weak 
solutions of the Euler equation, the connection between this limit and corresponding solutions of 
the Kelvin Hclmholtz equation remains a fully open problem. Moreover with the Prandlt-Munk 
vortex described in Section [5] one can show the existence of situations where such convergence does 
not hold. 

• Eventually one could consider a viscous perturbation of the water waves problem, i.e. a Navier- 
Stokcs equation with a free boundary. This has already been studied (cf. for instance [S3]). One 
may conjecture that in this case, with the Rayleigh- Taylor sign condition at the boundary, the 
convergence to the water waves solution should be true. 

Acknowledgments. Ce travail a beneficie d'une aide de l'Agence Nationale de la Recherche portant la 
reference ANR-08-BLAN-0301-01 
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